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Abstract 

Coupled ordinary differential equations are derived for the distant grav- 
itational interaction of a compact object of mass M and charge Q with an 
initially straight, infinitely long, cosmic string of tension /i <C 1/G = 1, when 
the relative velocities are very low compared to the speed of light c = 1. 
(An intermediate result of this derivation is that any localized force F(t) 
on the string that is confined to a single plane perpendicular to the initial 
string configuration gives the intersection of the string with this plane — the 
point where the force is applied — the velocity F(t)/(2/Li).) The coupled equa- 
tions are then used to calculate the critical impact parameter 6crit(^o) for 
marginal gravitational capture as a function of the incident velocity vq. For 
Wo ^ (1 — Q'^/M'^y/^fi'^/^, so that the string acts relatively stiffly, 
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For (1-(52/M2)V3^2/3 ^ -t;g ^ I-Q'^/AP, so that the string acts essentially 
as a test string that stays nearly straight, 
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imits the critical impact parameter is found numerically 
to fit a simple algebraic combination of these two formulas to better than 
99.5% accuracy. 
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1 Introduction 

The gravitational interactions between compact objects and long cosmic strings 
have been calculated in several recent papers [|l], ||, ||, ^ in the test-string approxi- 
mation, in which one assumes that the dimensionless string tension /i (using units 
in which Newton's gravitational constant is G = 1) is much smaller than any other 
parameter whose ratio with /x is relevant. Then the string obeys the Nambu-Goto 
equations of motion 0, H, |^ in the gravitational field of the object and has negligible 
back reaction on the motion of the object. 

It is indeed a good approximation for cosmic strings that /i <S 1, but I shall show 
that at low initial relative velocities vo and near the critical impact parameter, the 
ratio h/vq is important, so the test-string approximation is accurate only if the 
velocity is large in comparison to /i^/^. For velocities very low compared with this, 
the string acts relatively stiffly, so that the motion of the object responds more than 
the motion of the string. 

Here I shall derive coupled ordinary differential equations of motion for an ini- 
tially straight, infinitely long, cosmic string interacting gravitationally with a com- 
pact object of mass M and charge Q (i.e., one whose size is negligible in comparison 
with the distance between the object and the string), using the approximations 
yU <^ 1 and f <^ 1 but allowing yu/fo to be arbitrary. The low velocity allows 
the part of the string nearest the object to remain nearly straight, and then the 
force F on the string is a simple function of the string-object separation and relative 
velocity. On a much larger scale this force is nearly localized and hence gives the 
point on the string nearest the object the velocity F/(2/i), as we shall see. 

For fixed M, Q, and /z, after taking out Poincare transformations there is a 
two-parameter set of solutions that can be characterized by the initial speed Vq of 
the object relative to the string (when they are infinitely far away in the infinite 
past) and by the impact parameter b (how far the object would miss the string if 
both moved uniformly with no interactions). For each vq there is a critical impact 
parameter, 6crit, such that if 6 < bcrit the object becomes gravitationally captured and 
bound to the string, but if 6 > fecrit, the object and string scatter and asymptotically 
move freely of each other. In this paper the coupled equations of motion of the 
object and the string will be solved to give the function &crit('^o)) analytically when 
the ratio /i/f o is either very large or very small, and numerically otherwise. In the 
latter case simple algebraic combinations of the analytic results at the two extremes 
reproduce the numerical results over the entire range to within 1%, or even within 
0.43% if two arbitrary exponents are chosen appropriately. 



2 Gravitational force between an object and a 
straight string 

Strictly speaking, the gravitational interaction between an object and a string 
should be done in terms of curved spacetime. However, when fi <^ 1 and when one 
is in the weak-field regime of the object (which is where the string is, under the 
assumption that the object size is much smaller than its distance to the string), one 
can do an analysis in terms of fictitious forces acting in a fictitious fiat background 
spacetime, just as one does for Newtonian gravity (though one need not assume that 
the gravitational field of the string is Newtonian or that the object velocity is low 
to do this). 

We assume that the string is infinitely long and is nearly straight in the region 
where it is fairly near its minimum separation from the object, with no oscillations 
coming in along the string from infinity. Given the assumed initial straight, static 
configuration of the string, this assumption of approximate straightness of the part 
of the string nearest the object for all times is justified by the assumption of low 
transverse velocities, since the fact that the disturbances of the string propagate up 
and down it, away from their source, at the speed of light, means that the length 
of string that is disturbed is much greater than the transverse disturbance, so that 
the angular bending of the string is small, even if the linear bending becomes large. 

Then we can go to the frame in which the part of the string nearest the object is 
at rest along the 2;-axis and the object is in the x-y plane with velocity of magnitude 
Vr in the positive x-direction, the relative velocity of the object with respect to the 
piece of string nearest the object. To linear order in M and Q, the gravitational 
force of the string on the object is zero when Vr = 0. However, the string produces 
a conical spacetime with a deficit angle 27r(4/i) p, so when this is fiattened out in 
the fictitious fiat background spacetime with azimuthal angle ip that ranges from 
to 277, a moving object appears to be bending toward the string at the rate of A/j, 
times the rate at which ip is changing, i.e., at the rate of 4/i(y9 = AfiVry/r"^, where 



x'^ + y"^ (1) 

is the distance from the object to the string (ignoring relative corrections of order // 
that depend on how the conical spacetime is fiattened out to produce the fictitious 
fiat spacetime). 

When this bending rate is multiplied by the momentum p = M'jrVr of the object, 
where 7^ = l/Jl — v^, one gets the part of the transverse fictitious force on the 
object that is linear in M, A^M'jrv'^y/r'^. The longitudinal fictitious force (in the 
direction of the object velocity v^ in the instantaneous frame of the piece of string 
closest to the object, using boldface symbols to denote spatial vectors in the two- 
dimensional plane perpendicular to the string in the fictitious fiat spacetime) that 
is linear in M depends on the fiattening procedure used to generate the fictitious 



flat metric. One can choose it so that the total fictitious force on the object that is 
linear in M is 

Flinear = -^^JiM-frVlv / T^ . (2) 

By Newton's third law or conservation of momentum, there will be an equal and 
opposite fictitious force on the string to describe its motion in the fictitious flat 
spacetime. 

To explain and justify this fictitious force in greater detail, consider the geodesic 
motion of a test mass M in the conical spacetime of a cosmic string f^ , 

ds^ = -dt^ + dz^ + dp^ + (1 - ^.^fp^dif^ 

= -dt^ + dz^ + r-^^{dr^ + r^d^^) = -dt^ + dz^ + e^'^{dr^ + r^d^^), (3) 

where 

r= [(l-4/i)p]i/(i-^'^) (4) 

and 

= — 4/ilnr (5) 

(not to be confused with the azimuthal angle cp). 

One can go to the frame in which the object moves in the plane 2; = 0, in which 
case the motion is in an ultrastatic (2 + l)-dimensional spacetime with conformally 
fiat spatial sections, whose general metric form is 

ds^ = -dt^ + e^'*'5ijdx'dx^. (6) 

Timelike geodesies of such metrics have constant 

dt , ^ 

E = M— (7 

dr 

and constant 

P' = M^e^H,, — — = E^- M\ (8) 

dr dr 

One can also readily show that the spacetime geodesic equation implies that the 

spatial trajectory is itself a geodesic of the spatial metric e^'^5ijdx^dx^ and obeys the 

equation 

where a is the spatial distance along the trajectory, 

da"^ = 6ijdx'dx\ (10) 

in the fictitious fiat spacetime metric 

ds'^ = —dtl + Sijdx^dxK (11) 



One can reproduce the spatial geodesic of the confornially-flat spatial metric 
e^'^dijdx'^dx^ by a fictitious spatial force 

Y = ^ = M^,vlVci> + fwr (12) 

at I, 

in the fictitious fiat spacetime, where 

p = M-ir^r = M-irdr/dt\, (13) 

is the fictitious fiat spacetime momentum of the particle of mass M, of fiat spacetime 
position vector r (with components x* that are raised and lowered by the fictitious 
fiat spatial metric Sijdx'^dx^ in this fictitious force analysis), of spatial velocity v^ = 

dr/dt\, with magnitude Vr = y/Vr ■ v^ = \l 5ijVlvl, and of 7^ = l/y 1 — v"^. 

In the expression above for the fictitious force F, one can have an arbitrary 
function / multiplying the velocity v,, if one just wants to reproduce the spatial 
trajectory of the geodesic of the actual metric (H). If one also wants to reproduce 
the time dependence of the trajectory, so x*(tt,) under the fictitious force in the 
fictitious metric has the same dependence on t^ as x'^t) does on t for the geodesic 
in the actual metric, then one needs 

/ = -M7,(72 + l)v,-V0. (14) 

However, for the purposes of this paper, the precise time dependence of the 
trajectory of an object passing a cosmic string is not so important as the trajectory 
itself, so I shall set / = 0. Then if one uses the form of (f) given by Eq. (^ for the 
cosmic string, one gets the fictitious force given by Eq. (0). 

To quadratic order in M and Q, A. G. Smith |P calculated from the effect of the 
deficit angle on the three-dimensional Laplacian for the Newtonian gravitational 
and electrostatic fields of the object that the string exerts an attractive force (if 
M^ > Q"^) of magnitude 7r/i(M^ — Q^)/(4r^) on the object when the relative velocity 
Vr can be neglected. (When v^ is not negligible in comparison with unity, the 
quadratic force is overwhelmed by the linear force, so one does not need the high- 
velocity corrections to the quadratic force.) In two-dimensional vectorial form, 

Fquadratic = -T^ ^l{M^ - Q^ / {^T^) . (15) 

The conservation of momentum implies that the object exerts an equal but op- 
posite force on the string, and this indeed gives precisely the bending of a distant 
static string in the Reissner-Nordstrom metric ||T0|| . 

When we add the fictitious forces that are linear and quadratic in M and Q, we 
get a total fictitious force on the object that, in the limit of a small dimensionless 
string tension /i <C 1 and in the limit of a large separation r ^ M, is approximately 

Fm = -^ -^^ ■ (16) 



The first term is valid for all velocities v^ but neglects additional terms in the 
direction of v^, and the second term, which is only important if the velocity is very 
small, neglects corrections when Vr is not small. 

3 Motion of a long string with a localized force 

In the fictitious flat spacetime, the response of the compact object to the fictitious 
force is simply the relativistic version of Newton's second law, the first equality of 
Eq. ([T2|). However, the response of the string to the equal and opposite fictitious 
force needs to be worked out. 

As noted above, the low-transverse-velocity assumption means that the length 
of string that is bent (which develops at the speed of light) is much greater than the 
transverse bending (which develops at the transverse velocity of the string, much 
less than the speed of light by assumption). Hence, one can look at the bending on a 
scale much larger than that of the transverse separation between the string and the 
object and yet much smaller than that of the length of string that is bent. On this 
scale, the force on the string appears to be localized at the position on the string 
nearest the object and is given to good approximation by Eq. (|l^) in terms of the 
string-object transverse separation r and the relative velocity v^. 

Let us do a general analysis of the effect of such a localized force on an initially 
static, straight cosmic string, without, in the general analysis, assuming that the 
transverse velocity of the string is small, even though in our particular application 
we need this assumption to justify our assumption that the force is effectively local 
and is given by Eq. (|16]). 

Suppose a time- dependent horizontal force in the x-y plane is applied at z = 
to an infinite string initially lying at rest along the z-axis (vertically), with no other 
forces on it. This force causes equal horizontal disturbances to go up (for z > 0) and 
down (for z <Q) the string at the speed of light (propagating away from where they 
are generated by the time- dependent force, at z = 0; assuming that the infinite string 
was initially at rest and has no forces on it at z 7^ implies that no disturbances 
are coming in to z = 0): 

r,(t, z) = e{z)v^{t -z)+ 9{-z)r±{t + z), (17) 

where 6{z) is the Heaviside step function of z and each of the two occurrences of 
r^ represents the same function of the argument inside the following bracket (either 
t — z for the first term, which is nonzero for 2; > where 6{z) = 1 and 0{—z) = 0, or 
t + z for the second term, which is nonzero for z < where 0{—z) = 1 and 9{z) = 0). 
Focus on the infinitesimal piece of string just above z = a.t some arbitrary time 
t. Its tilt or slope away from vertical is 

drsit,z) drsit,z) . 



where the overdot in the last term denotes a derivative with respect to the unhsted 
argument t — z evaluated ai z = 0. Denote the magnitude of this slope by 



Vs = \r±\ = \/r±-r± = tcinilj, (19) 

where ip is the instantaneous angle of the tilt from the vertical. 

This piece of string can be considered to be moving straight upward at the speed 
of light (unity), the velocity at which the pattern moves for z > 0, but locally 
any component of the velocity parallel to the string has no effect. To calculate the 
energy and momentum carried by a piece of string, what is relevant is the transverse 
velocity, the component perpendicular to the string. Here the infinitesimal piece of 
string is tilted at an angle ip from the vertical, so the transverse velocity is tilted at 
an angle tt/2 — ip from the vertical direction along which the pattern speed is unity. 
Therefore, the magnitude of the transverse velocity is 

Vt = sinip = ^ , (20) 



^s 



and the associated relativistic 7 factor is 

lt= I — = sec V^ = y^l + vl- (21) 

The length of this infinitesimal piece of string is 

dl = {seci>)dz = ^l + vl dz, (22) 

and its energy is 

dE = H'jtdl = /i(sec^ ^)dz = /i(l + vl)dz. (23) 

The transverse momentum (at the angle ip from the horizontal x-y plane) is the 
transverse velocity times the energy, with magnitude 

dpt = VtdE = /i(sin-?/' sec^ il))dz = nVgJl + v1 dz. (24) 

The vertical component of this momentum is 

dpz = sintjjdpt = iJ.{ta,n'^ ip)dz = fxv'^dz, (25) 

and the horizontal component has magnitude 

dph = cos ipdpt = fj,(tanilj)dz = fiVgdz (26) 

and is in the direction of r^. 



This infinitesimal piece of the displaced-string pattern that is moving upward at 
the speed of hght from the force at 2; = takes a time dt = dz to be generated. 
A mirror-image pattern is moving downward at the speed of hght from the force 
at z = 0, and during the same infinitesimal time an infinitesimal piece of string 
with the same magnitude of dz is generated moving downward just below 2; = 0. 
The vertical components of these two pieces are opposite and hence cancel, but 
the horizontal components are in the same direction and so add to make the two 
infinitesimal pieces of string generated during the time dt have a total momentum 
of 

(ip = 2fir^dt = 2/i(ir^ = 2fidrs, (27) 

where here and henceforth I shall use r^ for rs{t, 0) at 2; = 0, the horizontal position 
of the string in the x-y plane where the force is applied. 

Since the string was assumed to start off from rest along the 2;-axis, at Vsit, z) = 
for large negative t, the total momentum of the string is simply 



and the force on the string is 



Ps = 2/xr„ (28) 



F. = ^ = 2/ir,. (29) 



The string will react back with an equal and opposite force on whatever is forcing 
it to move horizontally. It is interesting that this reaction force is precisely an 
idealized friction force directly proportional to velocity. The fact that momentum is 
carried away along the string at the speed of light, never to return (since the string 
is assumed to be infinite), is apparently what allows an initially straight, infinitely 
long string to act as a perfect source of friction. 

One application of Eq. (^) is to give a very simple derivation of the net displace- 
ment of an initially straight string by a gravitating object of mass M that passes by 
with large impact parameter b ^ M/vq, so that the fictitious force term linear in 
M, the first term of Eq. (ITBp, dominates over the force term quadratic in M and Q, 



the second term of Eq. ([TB|). Then the object trajectory bends by half the deficit 



angle, i.e., by an angle Atc/x <ti 1, so its momentum changes by 

Ap = ATTfip = Att^M'JoVo. (30) 

By Newton's third law, the string's momentum changes by the opposite amount, so 
Eq. (|28| ) implies that the string gets displaced by 

Ar = Ap/{2fi) = 27cM-^oVo, (31) 

just as De Villiers and Frolov found p|. 
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If one inverts Eq. ( P^ ) to get 

*• = !• '^2' 

one notices the curious fact that apparently a force larger than twice the string 
tension fi can move a string faster than the speed of light, seemingly a violation 
of special relativity, despite the fact that a completely (special) relativistic analysis 
was used in the derivation of the relationship between horizontal force and velocity. 
However, one can make two comments about this paradox: 

First, for all finite horizontal speeds Vg of the string, even 'superluminal' Vg > 1, 
the physical transverse speeds Vt of the string, given by Eq. ( ]20|) above, are less than 
the speed of light. The horizontal speed Vg of the kink in the string at z = where 
the force is applied is somewhat analogous to the speed of the point where the two 
blades of a pair of scissors intersect. Even though the two blades are constrained 
to move slower than the speed of light, there is nothing in special relativity that 
constrains the intersection point to move slower than light. 

One might object that the scissors are more nearly analogous to a pair of inter- 
secting strings both moving in the same plane, whose intersection point can indeed 
move faster than light (at least if one assumes no interactions at the intersection, so 
that the two strings each move freely) . In this case it seems to be the prior existence 
of each of the two strings on both sides of the intersection that allows the intersec- 
tion to move faster than light while neither string does so. However, for the single 
string being considered in this paper, whose kink at z = moves horizontally faster 
than light if the force is greater than twice the string tension, there is no pre-existing 
string beyond the kink, so it seems that something physical must be moving faster 
than light. 

The partial answer to this worry is that for a horizontal force greater than 2/i, 
string is being created faster than light, but it is not being moved faster than light 
(in the transverse direction, the only motion that has a physical meaning for a 
longitudinal-boost-invariant string of the type being considered here). This is some- 
what analogous to the possible nearly-instantaneous freezing of a lake, in which the 
water-ice boundary can in principle move faster than the speed of light as ice is 
created (but not moved) faster than light. 

Another analogy would be electron-positron pair creation in a hypothetical im- 
ploding cylindrical electromagnetic wave that builds up a huge longitudinal electric 
field, say between two capacitor plates at the ends of the cylinder. The pair creation 
can in principle discharge the plates much faster than the time taken for light to go 
from one plate to the other. If one draws the Feynman diagrams of the pair creation 
and annihilation that accomplish this process, with positrons pictured as electrons 
moving backward in time, one can get electron lines that zigzag forward and back- 
ward between pair creations and annihilations and have a net spatial motion faster 
than light between the two plates. One might call this a tachyonic lightning bolt. 
And yet there would be no violation of special relativity with signals traveling faster 
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than light in this hypothetical process, since it would be the imploding electromag- 
netic wave that locally causes the process, with no causes traveling faster than the 
speed of light. 

Considering the creation of string faster than the speed of light leads on to the 
second comment, which is that although a string with a kink having a horizontal ve- 
locity Vg > 1 does not seem to violate the principles of causality in special relativity 
(signals not traveling faster than light), it does seem to be associated with an in- 
stability analogous to tachyonic behavior. The example of the nearly-instantaneous 
freezing of a lake illustrates the instability of the liquid phase of water at low temper- 
atures, and the example of the pair creation in the imploding electromagnetic field 
illustrates the instability of the zero-current 'vacuum' of the Dirac electron-positron 
field in the presence of a strong electric field. Similarly, tachyonic negative mass- 
squared terms in relativistic wave equations, such as the Klein-Gordon equation, do 
not lead to violations of causality (since disturbances do not propagate faster than 
light whatever the sign of the mass-squared term) but instead lead to instabilities 
of exponential growth. 

In the case of applying a strong force to a string, this can lead to an instability 
in the production of strings. For example, suppose the force is that of a gauge field 
on a massless charged particle attached to the string at the kink aX z = 0. (If the 
particle were massive, its inertia would prevent the kink from moving faster than 
light, so I shall assume the contrary to continue the argument.) But then if the force 
were greater than /i, the energy that a pair of oppositely charged massless particles 
extract from the gauge field in being created and separated would be greater than 
the energy needed to create the string joining them. This would lead to the rapid 
pair production of pieces of string with these oppositely charged gauge particles at 
the ends in the strong gauge field, an instability that would presumably rapidly 
discharge the gauge field until the force on a charged particle would be no greater 
than /i. 

Therefore, if the force is that of a stable gauge field applied to a massless charged 
particle on the string, it appears that the maximum force is the string tension /x, 
and therefore the maximum horizontal velocity of the string kink where the particle 
is attached is one-half the speed of light (relative to the distant parts of the string 
that are still at rest along the z-axis according to the assumption that the entire 
string was originally at rest along the z-axis and that the force is applied only at 
z = 0). 
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4 Coupled motion of a slowly moving compact ob- 
ject and a nearly straight infinitely long string 

Now let us combine the expression (|16D for the force between a compact object 
and a distant nearly straight long string with Newton's second law for the response 
of the object and with Eq. (|32|) for the response of the string, when the relative 
velocities are all small compared with the speed of light, so that Eq. (|^) indeed 
applies for giving the force and so that the force is effectively localized (relative to 
the length of string that is bent) in order for Eq. (|3^) to apply as well. Use a 
fictitious flat spacetime metric with Cartesian coordinates chosen so that the string 
is initially at rest along the 2;-axis and the compact object moves in the x-y plane 
(z = 0), and use boldface letters to denote vectors in this plane. (The same letters 
when not bold will denote the magnitudes of the corresponding boldface vectors.) 
In particular, let tm denote the position of the compact object (with mass M and 
charge Q) relative to the origin (x, y) = (0, 0), let r^ denote the position of the string 
at z = relative to the origin, and let 



r — I'm — ^s 



(33) 



denote the position of the compact object relative to the string. 

Then the nonrelativistic coupled equations of motion for the compact object and 
the string become 



Mr 



M 



-2firs 



TTij{M^ - Q2)r 4/iM(r ■ r)r 



47^3 



It is now convenient to use the unit of time 

M 



T 



2/x 



and the unit of length 



J^ii. = 



16/i I M2 



1/3 



M, 



(34) 



(35) 



(36) 



which we henceforth define to be unity (unless otherwise specified). In terms of 
these units the speed of light is no longer unity but is 



vr 



y^ 



Q' 



-1/3 



-^ u 



2'' \~ M2^ 
Setting Tu = Lu = I, Eqs. (|3^ become 



TC 



y^- 



M2 



-1/3 



> 1. 



(37) 



TAf 



r 4yu(r ■ r)r 



(38) 
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At an arbitrary time tj, the freely specifiable initial conditions are rs(tj), rM{ti), 
and rM{ti) in the two-dimensional x-y plane, a total of six parameters. If one uses 
the Euclidean invariance of the x-y plane to take out translations and rotations, one 
is left with three Euclidean-invariant initial conditions, for example the object-string 
separation distance r, the object speed 



V = ItmI = \/rM -I'M, (39) 

and the angle ip between the string-object separation vector r and the object velocity 
vector 

V = tm, (40) 

defined so that 

r ■ V = rvcosip- (41) 

Instead oi ip, it is sometimes convenient below to use 

C = cosip = ^-^ (42) 

rv 

or 



^b 1 — C rv — Y-w ,, , 

X = sin- = A/ = W . 43 



Note that the v given by Eq. (|39D as the speed of the object in the frame in 
which the string was initially at rest, and used with this meaning henceforth in this 
paper, is not the same as the Vr used in Section 2, e.g., in Eq. ([T6|), which denoted 
|r|, the instantaneous relative speed of the object and the string, and it is also not 
the same as the Vg used in Section 3, which by Eq. ([T9|) denoted what we are now 
calling |rs|, the speed of the point on the string nearest the object, in the frame in 
which the string was initially at rest. 

If one takes out the time-translation invariance, one is left with only two initial 
conditions, which can be the initial speed wq at t = — oo and the impact parameter 
b = rsinip at t = — oo, where r is infinite and ip is — tt. There Eqs. (|38| ) imply 
that Ts = 0, and one can choose the initial string position to be r^ = (the string 
initially at rest at Xg = Vs = 0) and the object to be coming in from xm = — oo with 
yM = b and with velocity vo in the positive x-direction initially at t = — cxd. These 
initial conditions and the first of Eqs. (|38|) imply that 



r, = Vo - i-M = Vo - V. (44) 

Inserting this into the second of Eqs. (|38D gives a single second-order differential 
equation for the object vector position rM{t), although the second term in the right- 
most expression of Eqs. (^) makes this single equation nonlinear in tm- 

As the object comes in from infinity, it will be pulled toward the string, so that 
at first it will speed up and bend downward, and yM will decrease below its initial 
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value of b. Eventually the object will either be captured by the string (so that r and 
V tend to zero), or it scatter so that r increases to infinity again, depending on the 
initial velocity vq and the impact parameter b. For each vo, the boundary between 
these two qualitatively different asymptotic behaviors is given by the critical impact 
parameter 6crit(wo): for b < &crit(wo) the object is captured, and for b > &crit(^o) the 
object scatters to r = cxo with v asymptotically greater than zero. For b = bcnt{vo), 
the marginal capture case, r reaches a minimum and then slowly increases back to 
oo, but with V asymptotically decreasing toward zero, so that the object essentially 
has precisely the escape velocity from the string. 

The main goal of this paper, besides deriving the coupled equations of motion 
(|34D or (|38D given above, is to calculate the critical impact parameter fecritl'^o) as 
a function of the initial velocity vq when this velocity is low compared with the 
speed of light. I shall derive one approximate formula when L„/T„ -C fo <C c, the 
regime in which the string acts essentially as a test string at the critical impact 
parameter (at least until r reaches its minimum value; during the late stages when 
r slowly increases back to oo, the back reaction of the string on the object becomes 
important), and another approximate formula when vq ^ L^/T^ ^ c, the regime 
in which the string acts stiffly and is only slightly perturbed by the force from the 
object (though this small perturbation is essential for causing the object to lose 
enough energy to become bound). 

In this paper I shall focus on the case in which the ratio of the second term to 
the first term on the extreme right hand side of Eqs. (^8]), 4/ir(r ■ r), is much less 
than unity, so we shall omit the last term of Eqs. (^81). For the dominant part of the 
scattering, r is of the order of magnitude of the impact parameter b, and r-r is of the 
order of magnitude of the initial velocity Vq, so this ratio is of the order of magnitude 
of fibvQ (in units with Tu = Lu = 1). For vq <^ Lu/T^, we shall find below that 
focrit oc f , so for very small Wq, certainly /i6t>Q <^ 1 for the case of marginal capture. 
For Vq of order unity (Lu/Tu), the truncated Eqs. (|38D (with the last term omitted) 
have no small or large parameters, so we expect that then focrit is also of order unity 
(Lu), giving fibcritVo ~ /^ < 1- For vo > 1, fecrit ~ '^Vq'^^'^, so /ifecritfo ~ /^^o''^- This 
is also small compared to unity when vq -C fi'"^^^ ~ (1 — Q'^/M'^Y^^c. 

Therefore, for calculating the critical impact parameter 6crit('yo)) it appears to 
be a good approximation to neglect the last term of Eqs. ( p8| ) for velocities low 
compared with the speed of light times the cube root of (1 — Q^/M"^). Unless the 
object is significantly charged (e.g., nearly extremely charged if the object is a black 
hole), the limitation of considering only velocities low compared with the speed of 
light is sufficient for dropping the last term of Eqs. (^) when one calculates 6crit('?^o)- 
It is apparently not sufficient when {l — Q'^/M'^Y^^c <^ vq <^ c, but the calculation of 
bciit{vo) in this regime for nearly extremely charged objects will be left as an exercise 
for the reader. 

When we drop the last term of Eqs. (|38|) , use Eq. (|4^ to solve the first part of 
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these equations, and define 

R = tm - vo, (45) 

then what remains of Eqs. (|38D becomes the single second-order vector equation 



R = -^1±5-. (46) 

IR + RP ^ ' 



In terms of R, one has 



Tm = vo + R, (47) 

r, = vo - R, (48) 

r = R + R, (49) 

V = vm = R. (50) 

Alternatively, Eq. (|46|) can be written as two first-order vector equations. 






(51) 



^ = --r (52) 

If we instead go to the three Euclidean-invariant quantities r, w , and either C = 



cosip = r ■ v/(rf ) or X = sin {'ip/2) = w (1 — C)/2, one gets the three first-order 
scalar equations 

r = Cv--^, (53) 



and either 



or 



V = -^, (54) 

C=il- C') (^ - ^) (55) 

X = -Xil-X')(--^). (56) 

Another choice of variables is motivated by the fact that in the stiff-string regime, 
where r'^v ^ 1, the equations are approximately those of a unit-mass Keplerian 
particle orbiting a stationary string that exerts a unit inverse-square force on the 
particle, r ^ v and v = — r/r^. For such a particle, conserved quantities are the 
energy, 

E^\v'--. (57) 

2 r 



the angular momentum. 



L = rvsmi} = rv\/l - C^ = 2rvXy/l - X^, (58) 
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one-eighth the cube of the angular momentum, 



h = \L' = r\\X^-Xy\ (59) 

8 



the eccentricity of the orbit, 



e = Vl + 2EL2, (60) 

and one-fourth the excess of the square of the eccentricity over unity, 

g = i(e'-l) = -EL" = {r\^-2rv^)X\l-X^) = [{rv^ -1)^ -1]{X^ -X^). (61) 

Only two of these quantities are independent, and to get fairly simple differential 
equations below, I shall focus on g and h. In terms of them and X^ — X^ = 
(1/4) siia^ ip, one may readily solve algebraically for r and v. 

/l2/3 

(62) 



VX2 - X4(VX2 - X4 + VX2 -X^ + g ) ' 



V = /l-V3(VX2 - X4 + ^X2 - X4 + ^ ). (63) 

Because the string is in actuality not completely stiff but moves toward the 
object, the quantities E, L, h, e, and g are not precisely conserved in the coupled 
object-string equations but instead are all decreasing with time, at the rates 

E = -^, (64) 

L = -^, (65) 

h = -^ = -3h-\X^ - X^)3/2(VX2-X4 + ^X^-X^ + g f, (66) 

e = -—^{rv'-l) = -^C, (67) 

er* e 

[2 1 

g = -(rv^-l) = — v^C = 2v^XX 



= -2/i-2(x2 - X^)=^/YX2 -X^ + g (VX2 - X4 + y'x2 - X^ + ^ )^ (68) 
Now we can insert Eq. ( |B^ and Eq. ( |BBD into Eq. (^B|) to get 

(X2 - X4)3/2 



X = -^ ^^ ^ V^X2-X4 + ^(VX2-X4 + y^X2-X4 + ^). (69) 
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Then dividing g and hhy X gives two coupled first-order differential equations for 
the coupled motion of the object and the string: 

dg 2 



^^ ^X{^X^-X^ + ^X^-X' + g)\ (70) 



dh_ ^ 3X(VX2 - X4 + yx^ -X^ + gf 

dX ~ ^X'^ -X^ + g ' ^ ^ 

At t = — oo, one has the object approaching the string, and so the angle V 

between the object velocity v and the infinitely-long string-object separation vector 

r is vr, giving X = sin {1^/2) = 1 initially. The energy is then Eq = Vq/2, and the 

angular momentum is then Lq = bvo, so 

9iX = 1) = IeoLI = Ife^o', (72) 

h{X = 1) = Ul = Ifo^ol (73) 

Then as X decreases, so do g and h. If the object scatters and remains unbound, 
the energy and hence g remain positive as the object returns to infinitely large r, 
but now with the object moving away from the string, so v and r asymptotically 
approach the same direction, giving ip = and hence X = 0. If the object becomes 
bound to the string, E = v'^/2 — 1/r becomes negative, and hence also g. Eq. ([7l|) 
becomes singular where a/X^ — X^ + g goes to zero, but one can switch back to 



Eqs. (|53D-(|55D or (0) to continue the calculation, as X and a/X^ — X"' -|- g reverse 
sign and the evolution continues, with X oscillating as the object spirals in to join 
the string after a finite time. 

In the intermediate case in which the impact parameter b has the critical value 
&crit('i^o); E and g tend to zero as the object slowly moves to infinite separation r 
with essentially just the escape velocity. In this case, like the scattering case in 
which g stays positive, the object velocity v and the string-object separation vector 
r asymptotically become parallel, giving X = at t = cxd. 

It is now easy to calculate a one-parameter family of vq's and their corresponding 
bcTitivoYs using Eqs. (^) and (0), without requiring a shooting method from initial 
conditions at t = — oo or X = 1. Instead, one evolves backward in time from t = -|-cx) 
or X = 0, where one sets g{X = 0) = but uses 

rj = h{X = 0) = ll^t = oo) (74) 

8 

as the single free parameter for determining the one-parameter family of critical 

solutions. Each solution is evolved to X = 1, where one evaluates g{ri,X = 1) and 

h{ri,X = 1) and then inverts Eqs. (|72D and (^) to get 

7;o(r7) = g{v, X = l)'/'h{r], X = l)-l/^ (75) 

bcritivoiv)) = 2g{v,X = l)-'/'h{7],X = l)^/^ (76) 
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5 The critical impact parameter in the stiff-string 
regime, ^o < Lu/Tu 



Now we shall find approximate analytic solutions of the evolution Eqs. ( ffOD and 
( ffTD in the marginally bound case {g{X = 0) = or E{t = oo) = 0) in the two 
limiting cases in which rj = h{X = 0) = L^{t = oo)/8 is either very large or very 
small. 

For 1] ^ 1, the factor of h in the denominator of Eq. ([70|) means that to lowest 
order g{X) stays very small and can be neglected on the right hand sides of Eqs. 
(|70D and (|Tl|). Then one gets 



9{X) 



64:7] 



(12'0-8sin2^ + sin4^) 



— [3 sin-i X + (-3 - 2X2 ^ 24X^ - 16X^)VX^ - X% 

oTj 



h{X) ^r7 + -(2zA-sin2z^) = rj + -[sm~^ X + {-1 + 2X^)VX^~^X^]. 
8 2 



This is sufficiently accurate to give 

h{v,X = l) 



4:7] 7f 



(77) 
(78) 

(79) 



However, to get g{7], X = 1) with a relative error of only 0{ri~^) rather than 0{r]~^), 
one needs to insert the approximate expressions (77) and (|78|) into the right hand 
side of Eq. (|70D and evaluate it through the next-to-lowest order in 7]~^. This gives 



9{V,X 



37r 
IQ7] 



1 + - + 0(1) 



(80) 



If we insert these expressions into Eqs. ( [75D and (|76D to get fo(^) and 6crit('7) 
and then eliminate 77 to get hcntiyo), we find that 



5 



(81) 



Large 7] gives small v^ ^ \/3tt/{47]^/^), so Eq. ( PTf ) applies for vq -C Lu/Tu = 1. 
At the critical impact parameter, this is the stiff-string regime, where the string 
responds much less than the object to the force between them. If one restores the 
units Tu and Lu given by Eqs. (^) and (|36|), one obtains the formula given in the 
Abstract above for the stiff-string regime. 

One can give an order-of-magnitude derivation of the critical impact parameter to 

—7/5 

show that it indeed goes as Vq for wq -C 1. What happens is that the object comes 
in on a slightly hyperbolic, but nearly parabolic, orbit with initial energy Eq = Vq/2 
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(83) 



just slightly positive. Then as the object swings around the string, the small motion 
of the string toward the object drains just enough energy from the object that it 
goes back out on an asymptotically parabolic orbit with E approaching zero. 

In this analysis, use ~ for approximate numerical equality under the limits given 
(i.e., a; -C 1 implies that 1 + x ~ 1), and use ~ for order-of- magnitude equality that 
has the correct powers of small quantities (e.g., 2a; ~ x). 

With impact parameter b ^ 1/fo, the angular momentum L stays nearly con- 
stant at its initial value Lq = Bvq ^ 1 (as I shall show momentarily), so at the 
perispaggon, the point on the object orbit nearest the string, the minimal separa- 
tion Tm and the object speed Vm are related by r^fm ^ L ^ Lq = Bvq. Assuming 
that E = v'^/2 — 1/r <^ 1/6 (which for E r^ Eq requires h <^ l/fg, so we need 
1/fo <^ 6 ^ 1/fo and hence fo <^ 1), we get v'^/2 ^ l/r^ ^ vl/2. Combining this 
with TrnVm ~ Lq then gives 

Tm ^ \lI (82) 

2 

The time spent near perispaggon is tm ~ fmlvm ~ -^o- During this time, the an- 
gular momentum decreases at the rate given by Eqs. (^), L = —L/r^ ~ —Lq/t^ ~ 
— 1/Lq, and thus it decreases by a total amount roughly tm as large, or AL ~ — 1/Lq. 
The magnitude of this decrease is much less than Lq for Lq ^ 1 or 6 ^ l/'^O) as 
claimed above, so in this regime the angular momentum is approximately constant. 

The energy decreases at the rate given by Eqs. (0), E = — 1/r^ ~ — 1/r^ ~ 
1/Lq, and multiplying this by the time tm gives a total decrease —AE ~ ^/Lq ^ 
l/(6fo)^. For the final energy to be zero, this decrease needs to equal the initial 
energy, giving t>Q ~ i?o ~ —AE ~ l/(6t>o)^ or b = fecrit ~ "^o • For vq <^ 1, this is 
indeed within the allowed range 1/fo ^ &crit ^ Vq <^ l/fg of the approximations 
used above. 

One can easily get the correct coefficient (Gtt)^/^ of the Vq term for fecrit by 
calculating the time integral of the energy loss rate E = —1/r^ for a parabolic orbit 
of angular momentum Lq = Bvq, but that confirmation of the precise form of the 
first term of the right hand side of Eq. (^Tj) will be left as an exercise for the reader. 
Similarly, one should be able to get the second term on the right hand side of Eq. 
( |8lD from first-order perturbations due to the motion of the string and the resultant 
loss of the energy and angular momentum of the object. 

6 The critical impact parameter in the test-string 
regime, ^o > L^/Tu 

Now I shall consider the opposite case, rj = h{X = 0) = L^(t = cx))/8 <^ 1, 
which leads to vq ^ Lu/Tu = 1 and fecrit <^ L^ = I. This is the test-string limit, for 



which the lowest-order approximation is to ignore the string tension and the effect 
of the string on the object. The critical impact parameter in this limit was found in 
[|] to be 6crit — 2t'o after converting to the units used in this paper. Here I shall 
rederive this result and give the lowest-order correction in the small quantity 1/vq. 

In this test-string regime fo ^ 1, the qualitative picture is that initially the 
object moves with very nearly constant velocity, and the string gets bent toward 
the object as described in 0. In that paper, where the back reaction of the string 
on the object was totally ignored, the final situation, when the impact parameter h 
has the critical value focrit for marginal capture, is for the string to develop a kink 
that asymptotically trails the object at a distance of f q = &crit/2 behind it, an 
unstable equilibrium point where the force of the object on the string causes the 
string to follow the object at precisely the same speed as the object, by Eq. (|32D. 
In the complete test-string approximation, the deceleration of the object is ignored, 
so it continues moving forever, pulling an ever longer kink of string behind it. 

When one does include the back reaction effect of the string tension on the object 
but stays in the test-string regime oi vq ^ 1, the picture during the first stage of 
the evolution in the marginal capture case is very similar to that of the complete 
test-string approximation used in 0, so the position on the string nearest the object 
approaches quite near to what would be the unstable equilibrium separation behind 
the object if the object continued moving with constant velocity. But then in actual- 
ity the object slows down as the string weakly pulls back on it (somewhat analogous 
to the way a fish slows down when attached to a weak fishing line [|l3]). If it had 
been at the unstable equilibrium separation for the velocity of the object (i.e., the 
separation that would be the equilibrium one if the object continued moving with 
constant velocity), then the slowing down of the object would eventually lead to its 
falling into the string. For the object to be just marginally not falling into the string, 
the separation must be greater, by a precise tiny amount, than what would be the 
equilibrium separation at constant velocity. At this slightly greater separation, the 
force on the kink of the string trailing the object would be just slightly less than 
needed to move it at the speed of the object, so the string moves slightly slower 
than the object and hence trails behind it at an indefinitely increasing separation, 
gradually slowing down more and more. 

If the object had had an impact parameter h slightly greater than 6crit('yo); it 
would maintain a velocity bounded away from zero as its separation from the string 
increased indefinitely. On the other hand, if the object had b slightly smaller than 
&crit(f^o)) the object would eventually slow down more than the string, so that the 
separation would stop increasing and would henceforth decrease to zero as the object 
falls into the string. But in the marginal capture case with impact parameter h = 
^crit('yo); the separation r between the string and the object, although continuing to 
grow indefinitely, grows slowly enough that the object, as well as the string, continues 
to slow down just fast enough that it asymptotically approaches, but never reaches. 
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zero velocity. 

We can put this discussion into a more quantitative form in the following way: 
For 5f < X^ < 1, Eqs. ( 



I) and ([71]) become 

dg__ IQX^ + l2X^g + 0{X^) + 0{g^) 
dX~ 



dh_ 
dX 



h 



12X^ + 0(X') + 0(g^), 



with the critical solutions {g = and /i = r7>0atX = 0) 

h{X) =r] + AX^ + O(X^) + 0{7]'^). 



O 



X 



11' 



(84) 
(85) 

(86) 
(87) 



To get the time dependence of these solutions, we can use Eq. (|69|), which for 
g <C X^ <ti 1 becomes 



X 



4X^ + 3X'^g + 0{X^) + Oig"^ 
2h 



Inserting g{X) and h{X) from Eqs. ( p6D and (|87|) into this and integrating gives 

"3\ 



6X3 5 I ^ 



^9) 



dropping the irrelevant constant of integration. For large t this can be inverted to 
give 



X=(^ 



1/3 



In^t' 



45t r] \ t^ 



Then from Eqs. (H), (JS^), (H), and (H), one may get 

,(t)^f,^/3(6t)-/3(l + lln^), 
A 3 8 6t\ 

Alternatively, one can divide Eq. (|53D by Eq. (|^) to get 

-— = sec ^jJ — r^f ^ 1 — r'^v 
dv 



(90) 

(91) 

(92) 
(93) 
(94) 

(95) 
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for X = sin if)/ 2 <^ 1, where one may recall that ip is the angle between the velocity of 
the object and the string-object separation, so ip goes from vr initially (as the object 
heads toward the string) to finally for marginal capture (as the object heads away 
from the string with slower and slower velocity). When X can be neglected, Eq. 
(|95D leads to the series solution 

1 ~ 1^2 _ 1 5 ^ ^8 _ J_^ll ^ ^^14 Q( 17) (gg) 

r 2 20 160 8800 9856 ^ ^ ^ ^ > 

or 

v{t) ^ V^r-V2 ^ 1-2 ^ V|^_7/2 ^ 24 ^_5 ^ ^(^-13/2)^ (97) 

5 25 1375 

All of the series expansions given so far in this section are good only when t ^ 1, 
r ^ 1, and f <^ 1, in addition to X <^ 1 (i.e., when the object is moving very slowly 
nearly directly away from the string at very late times). When r^ <^ 1, as one 
integrates backward in time from the final condition (7 = and /i = r/ at X = 
(where t = oo, r = oo, and f = 0), the approximate second equality of Eq. ( P5| ) 
remains valid so long as X <^ 1 and gives, when v becomes very large and r becomes 
very small, 

i ^ ^1/2 _ 1^-1 _ l,r5/2 _ 15 4 _ 112^^-11/2 Q( -7) (gg) 

r 4 32 64 2048 ^ v ;, v ; 

up to terms that decrease roughly as the exponential of — (4/3)f^/^, so long as X 
remains negligible. Inverting this gives 

^(0 ^l-2 + \r + y + y^ + 0{r'') (99) 

for small r, again neglecting all contributions from nonzero X and terms exponential 
in -4/(3r3). 

To get the X-dependence of the critical solution when X <^ 1 and to prepare 
the way for the corrections when X is no longer small, divide Eq. (^) by Eq. (|5^) 
to get 

with the approximate equality applying when X ^ 1. 

When w <C 1 so that Eq. (p6| ) is applicable, Eq. ( |100[ ) becomes 



dX X 

dv V 



(101) 



which, when matched to the boundary condition that L = 1rvX\/\ — X^ is equal 

to 277-*^/^ at X = 0, gives 

1 , r?^/3 

X ^ ir^V3, ^ ^ (102) 

2 \/ It 
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for f ^ 1, r ^ 1, or X <^ rj^^^. 

At X ~ 1]^^^ <^ 1, f surpasses unity and thereafter becomes large, and r becomes 

small (until X gets near 1 and r becomes large again in this evolution backward in 

time). When f 3> 1 so that Eq. (^) is applicable, with rf^ ^ f^/^ ^ 1, Eq. (|100|) 

becomes 

dX X{l-X^)rv 

where the first approximate equality applies for all X ^ r]^/^ and the second ap- 
proximate equality applies for 7]^^^ ^ X ^ 1. 



Eq. ( |103[ ), along with Eq. (PHf), implies that once v gets large, X grows roughly 
exponentially with (2/3)f^/^. After X passes t;"^/^, dv/dX becomes small, and the 
large v changes very little during the rest of the backward evolution to the initial 
conditions at X = (t = — oo, v = Vq). Since this exponential growth of X with v 
starts at X ~ r]^^^ where f ~ 1 and ends at X ~ 1 where v ~ f o, one gets 

t;o~(-^lnr^)2/3 (104) 

for very small rj. This very slow growth of Vq with l/rj implies that when one solves 
Eqs. (^) and (^) numerically with the numerical initial conditions g{X = 0) = 
and h{X = 0) = r/ (final conditions in time for marginal capture), one must choose 
7] to be extremely tiny to get a large value of vo, and for the approximations of 
this section, one must assume by r^ ^ 1 not just that rj is more than an order 
of magnitude (say) smaller than unity, but that — In r^ is more than an order of 
magnitude (say) larger than unity. 

Once the marginal capture solution for r/ ^ 1 has been carried past X ~ 
(— Inr/)"^/^ where dv/dX becomes small, v stays nearly constant and one is in 
the test-string regime treated in 0, where one can to first approximation neglect 
the back reaction of the string on the object (the object's acceleration toward the 
string). In this paper, however, we want to treat this back reaction perturbatively 
in this test-string regime and get the first-order correction to the test-string critical 
impact parameter formula fecrit — 2fo given in |^. 

For this purpose it helps to change the dependent variables of integration from 
the g{X) and h{X) used in Eqs. (^) and ([7I| ) to two other variables that would 
be constant if the test-string approximation were completely accurate. One such 
variable is v, the magnitude of the velocity of the object (in the frame in which the 
string was initially at rest along the 2;-axis). A second such variable is motivated 
by the test-string Eq. (67) of |0, where the polar coordinate angle 9 defined just 
before Eq. (54) of that paper is essentially the same as ip here, where r there is 
the same as r here, where z given there just before Eq. (54) as r sin 6 is hence the 
same as r sin tp = 2rXy/l — X^ in the notation of the present paper, and where the 
constant L given by Eq. (50) (not to be confused with the angular momentum L 
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of the present paper) is the same as Vq in our present units Tu = 1 and Lu = I- 
Then Eq. (67) of B, plus the equation 6crit — 2L three hues below, imply that 



bi,,c:.Ar'X\l-X' 



Av^'X' 



(105) 



is a constant of motion in the test-string approximation, and furthermore that it is 
zero for marginal capture. If we use the fact that f ~ f g to replace Vq by v and then 
multiply the expression on the right hand side by v/4, we get the quantity 



w = r\{X^ - X^) - X' 



Av 



-X' 



h 



VX2 - X4 + VX^ -X^ + g 



Xl (106) 



Eqs. (70) and ([7T|) now imply that the alternative dependent variables v and w 
obey the differential equations 



dv 
dX 



X(l - 2X^ 



v{X^ - X^){X^ + w) 



1- 



X2-X4 



\ t;3(X2 + 



w) 



(107) 



dw 
dX 



X[X^-{l-2X^)w] 
^v^{X^-X^){X^ + w) 



X2-X4 
M^3(x2 + u;) 



(108) 



We can see from Eq. 
the way from X = 0, since sufficiently near there that h 



that w is not a convenient variable to integrate all 

Tj and g <^ 1, the last 
expression of Eqs. ( |106| ) implies that we have w ~ h/{2X), which diverges at X = 0. 
Alternatively, in the second-to-last expression of Eq. ( |106| ), one gets a divergence 
from the fact that L goes to the nonzero constant 2r]^^^ at X = 0, but v goes to zero 
there. However, from Eq. (p8|) for t]^^^ ^ X ^ 1, where f ^ 1 and r ^ t>~^/2 <^ i^ 
one can see that then w ~ {l/2)v~^^'^X'^, which is negligibly small for small enough 
X. E.g., for X ~ ?7^/^ ~ e~(2/3)i)o ^ where v ^ \ and r ~ 1, w ~ e"^^/^^*'o is 
exponentially small for large vq. Therefore, with only an exponentially small error, 
we can set w = and f :^ 1 as initial conditions in integrating Eqs. ( |107|) and ( |108|) 
from some very small X, such as r/^'^, which for practical purposes we can replace 
with with only an additional exponentially small error. 

Because Eqs. ( [L07| ) and (|108| ) have f ^/^ 2> 1 and v^l"^ ^ 1, respectively, in the 
denominators of the right hand sides, whereas the other factors are of order unity 
(except for harmless square-root singularities in 1 — X^ at X = 1), f and w change 
by only small amounts as one integrates to X = 1 to get the temporally initial 
conditions for marginal capture. To lowest order in 1/fo, one can thus neglect w on 
the right hand sides and rewrite Eqs. ( |107| ) and ( |108| ) as 



dv 



2X2 



2X2 



^^ Jt;(X2-X4) Jt;o(X2-X4) 



(109) 
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dw X^ X^ 



dX ^v\X^ - X4) ^vl{X^ - X^) 
Then Eq. (|110|) integrates to give 



(110) 



w{X = l) = ^^-l = \bl,,v, - 1 = V,''' + 0{vf). (Ill) 

From this we can readily solve for the critical impact parameter, in the units of 
Tu and of L„ given by Eqs. ( P3| ) and (PU]), for large Vq ^ Lu/Tu = 1: 

6erit(t^o) = 2vo'^^ + t;o-2 + 0(i;o"'/'). (112) 

When the units are restored, one gets the result given in the Abstract. 

— 1/2 

An order-of-magnitude argument for the Vq power law of the leading term 
was given in |^. Summarizing it very briefly in the units used here, if the object is 
incident with velocity f o 3> 1 and impact parameter b such that it is not captured, 
the distance of closest approach will be r ~ 6, the velocity at that point will be 
f ~ fo, the time when the string is roughly that close will be ~ r/v ~ b/vo, 
the length of string that bends significantly will be the speed of light times this 
or ~ cb/vo, the bending angle will be ~ l/(cr^) ~ l/(c6^), and hence the total 
transverse bending of the string will be the length of string significantly bent times 
this bending angle, or ~ l/{bvo)- Marginal capture requires that this total transverse 
bending be of the same order as the impact parameter, giving b = bait ~ Vv^- 

7 The critical impact parameter over the entire 
range of nonrelativistic velocities 

Now that we have found analytically Eq. (^) for the critical impact parameter 
at very low initial velocities, fo <^ Lu/Tu <^ c, and Eq. ( |112[ ) for bcntiyo) at high (but 



still nonrelativistic) velocities, L^jTu <^ f o <^ c, we would like to find it numerically 
for all nonrelativistic velocities, and ideally find a fairly simple formula that gives a 
good fit to 6crit('yo) for all v^. 

It was straightforward to program Mathematica 2.0 on my ancient NeXT to 
integrate Eqs. ( [70|) and (^) numerically from the temporally final conditions f? = 
and /i = r/ at X = back to the temporally initial conditions at X = and then 



use Eqs. (|75D and (^) to solve for t>o and &crit(^o)- It is slightly awkward that by 
this method one cannot straightforwardly choose the values of Vq for which one gets 
bcn\.{yiQ)^ but this is a small price to pay for not having to use a shooting method to 
find the critical or marginally bound solutions. (One could use a shooting method to 
get any particular f q desired, but for determining the characteristics of the function 
&crit(wo)5 this is hardly necessary.) 
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Another slight awkwardness is that to get large values of the initial velocity fo, 
one must use exponentially small values of t], and this can lead to some numerical 
inaccuracy, but this was not a serious problem for at least moderately large values 
of fo- 

I calculated over 50 pairs of (fo, ^crit('?^o))? ranging from (0.00104198,26910) 
through (for example) (0.834385,3.48671) to (357.89,0.105727). Larger values of 
vq require such extraordinarily small values of rj that it would be somewhat difficult 
by my numerical method to obtain them, but what I could easily attain extends 
well into the domain of validity of both Eq. (|8lD and Eq. (|112|) . 

The results from the numerical integrations fit extremely well to the truncated 
series expansions of Eq. ( pTD at very small vq and of Eq. ( |1 1 2| ) at large f o, confirming 
both the two coefficients and the two exponents of each of the two expressions (eight 
parameters calculated analytically and confirmed numerically). 

One can combine Eqs. (^) and (|112|) into the interpolation formula 



bcTit{vo) ~ B(^a,f3){vo) 
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(113) 



which depends not only on the eight parameters (four coefficients and four expo- 
nents) calculated analytically and given in Eqs. 



8T|) and 



but also on two new 
parameters, a and f5. Let us see how well this formula works for suitable a and (3. 

In order that the first two terms of the right hand side of Eq. ( pj) be the 
dominant terms of the right hand side of Eqs. ( |113| ) at fg ^ 1, one needs a > 4/3 
and /3 > 0, and in order that the first two terms of the right hand side of Eq. ( |112| ) 
be the dominant terms of the right hand side of Eqs. (|113|) at fo ^ 1, one needs 
a > 5/3 and /? > 0. The combination of these two pairs of inequalities gives what 
might be called Criterion A, a > 5/3 and (3 > 0. If one wants to avoid any deviations 
larger than the 0{vq) term of Eq. (^) at small Vq, one needs a > 8/3 and (3 > 2/3, 
and if one wants to avoid any deviations larger than the 0{vq ) term of Eq. (|112| ) 
at large fo, one needs a > 10/3 and (3 > 5/6. The combination of these latter two 
pairs of inequalities gives what might be called Criterion B, a > 10/3 and f3 > 5/6. 

If one wants a simple choice of the pair of fitting parameters (a, P) (e.g., a pair of 
integers). Criterion A motivated my first choice, (a,/3) = (2, 1). As a crude fit, this 
is not too bad, giving a maximum relative error of slightly over 10%, e.g. 10.382% 
error at (fo, &crit(fo)) = (0.940172,3.10142). My second choice, {a, (3) = (3, 1), is the 
smallest pair of integers that (almost) satisfies Criterion B. This gives 
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This approximate formula gives a much better bit, with a maximum relative er- 
ror of only about 0.85%, near (fo, &crit(fo)) = (1.79814,1.78528). (The position of 
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the maximuin relative error was not calculated to the six-digit figures given for 
(■^0) bcritivo)) here, but they are given to six digits so that one can see how accurate 
this, or some other, approximation is at this sample point, a point where the relative 
error was about as large as I could find given the six-digit apparent accuracy of the 
Mathematica calculation.) 

I then used the FindMinimum function of Mathematica to minimize, as a func- 
tion of the pair {a, j3), the sum of the 100th powers of the relative errors of -B(o,/3) (f o) 
over the first 45 pairs of values of (vq, ^critl'^o)) that I had calculated. This sum is 
an analytic function of (a,/?) whose minimum should be at a value of {a, (3) which 
is near the point where the maximum absolute value of the relative error for the 45 
points is minimized, since the sum of the 100th power of the relatives errors should 
be dominated by the largest relative error. 

Mathematica gave the minimization as occurring at {a, (3) = (4.00448, 1.36884), 
though the result would no doubt have been slightly different if I had used a different 
set of (fo, &crit('yo)) and a different power of the errors. Ideally one should use an 
infinite set of points and an infinite power of the errors, but I considered 45 points 
and a power of 100 as being high enough to give a good feel for the true minimum, 
over [a,f3), of the maximum of the magnitude of the relative error. The largest 
relative error I found for this value of (a,/9) was 0.427328%. 

By comparison, minimizing the sum of the squares (2nd powers) of the relative 
errors led to {a, (3) = (3.98412, 1.37903), which is not that much different even for 
this relatively small value of the power, though the maximum error I found for this 
choice of {a, (3), 0.47019%, is about 10% higher. 

A fairly simple rational approximation for the (a,/3) that minimizes the maxi- 
mum relative error of B(^a,/3){vo) is (4, 15/11), which then leads to 
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Here the number 4 in the exponents of the first term corresponds to a, and the num- 
ber 11 in the exponents of the second term corresponds to 15//?, so if one wanted 
different values of [a, (3), one could simply change these two numbers in this ex- 
pression. This approximation gives a maximum relative error of just under 0.43% 
(I found a maximum of 0.429991%), near {vo,bcrit{vo)) = (0.820853,3.54436). Thus 
this -5(4,15/11) (vo) is over 99% as good as -B(4.oo448, i. 36884) ('i'o) (has a maximum relative 
error less than 1.01 times the maximum relative error of -B(4.oo448,i.36884)('yo)) and is 
roughly twice as good an approximation (i.e., with roughly half the maximum rela- 
tive error) as the simpler i?(3^i)(t>o) given by Eq. ( |114|) that I guessed on my second 
attempt without doing any numerical error minimization. To get a relative accuracy 
better than roughly one part in 234 for 6crit('yo)5 one would need a different func- 
tional form than B(^a,f3){vo) given by Eq. (|113|) , but I won't pursue such refinements 
here. 
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—7/5 

Since the critical impact parameter 6crit goes as Vq for small velocity vq and 

— 1/2 

as Vq for large fo, one can see that if one multiplies focrit by fo to get the critical 
initial angular momentum Lcrit(fo) (per unit object mass M, and in units of Ll/T^ 
after dividing out by this M), this rises with both low and high vq and hence must 
have some minimum in between. By calculating several points near the minimum, 
I found that the minimum value is 

Linin ^ 2.9026, (116) 

which occurs near the sample point (fo, &crit(f^o)) = (0.834385,3.48671) given above. 
When one restores the units and multiplies by the object mass M to get the true 
minimum angular momentum of the object about the initial static, straight string 
configuration, one gets 



L^i, ^ 2.9026 



32^1 M2 



M^ ^ 1.9611/1^/3(1 - QV^')'^^^'- (117) 



A compact gravitating (and possibly electrically charged) object approaching from 
infinity toward an initially static, straight, infinitely long cosmic string with angular 
momentum (about the initial string axis) less than this would inevitably be captured 
by the string. 

8 Crude formula for the critical impact parame- 
ter for all velocities 

In the test-string approximation ||^, I calculated the critical impact parameter for 
velocities much lower than the speed of light and then a correction term linear in the 
velocity divided by the speed of light, getting Eq. (88) of that paper. In the case of 
a Schwarzschild black hole (Q = 0), I arbitrarily added a (relatively small) constant 
term to make the resulting formula agree with the known critical impact parameter 
at the speed of light ||l||, which is the same as the critical impact parameter for a 
massless point particle, ^JTJM. The resulting formula, Eq. (89) of 0, though just 
a guess by no means rigorously derived, is apparently a very good approximation to 
the actual critical impact parameter for a test string at all calculated velocities (i.e., 
up to a relativistic gamma- factor of 70 ~ 10) [^ 0, though it misses an infinite 
series of discontinuities that I predicted [|| occur at large values of 70 (none seen yet 
in the analyses of 0, Q up to 70 ~ 10). At velocities at all comparable to the speed 
of light, the critical impact parameter apparently must be calculated numerically by 
solving partial differential equations 0, ^ , which is a much more difficult numerical 
problem than using the ordinary differential equations of |^ and of the present paper 
that are only accurate at velocities very low compared with the speed of light. 



27 



Here I would like to give an analogue of Eq. (89) of [Q that incorporates (a) Eq. 
(|115|) as an excellent approximation to my numerical calculations at all velocities 
very low compared with the speed of light, (b) Eq. (88) of f^ for the term linear 
in velocity that is independent of the string tension, and also (c) the right limit for 
the critical impact parameter at the speed of light for a Reissner-Nordstrom black 
hole (arbitrary Q^ < M^). This last value is the minimum value of J —gee/ Qm in 
the Reissner-Nordstrom metric, which is 
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The resulting crude formula for all velocities for the critical impact parameter 
for the capture of a Reissner-Nordstrom black hole of mass M and charge Q by an 
initially straight, static cosmic string with dimensionless tension /i = G^ ^ 1 is 
that &crit('yo) is roughly the same as &guess('yo)) with 
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The first two terms dominate for f ^ 1 (where we have now returned to units in 
which the speed of light is c = 1). For this highly nonrelativistic velocity these two 
terms give essentially the same as Eq. ( |115| ), which I found numerically is accurate 
to greater than 99.5% accuracy for all nonrelativistic velocities when /x -C 1. The 
terms which were added to Eq. ( p,15| ) to get these two terms are essentially constants 
(up to correction factors that go to unity when /x is taken to zero) at Vq ^ jj?/^ and 
were designed to make these two terms vanish at fo = 1. The third term on the 
right hand side of Eq. ( |119| ) is the contribution of the term linear in vq calculated 
for the test string in [^ (though with the recognition that there might be other such 
terms missed in that calculation, and perhaps other terms that are independent of 
vq and that go as Vq ), corrected with a constant term to make this third term also 
vanish at vq = 1. The fourth (and last) term on the right hand side of Eq. ( p.l9| ) is 
independent of Vq and gives the critical impact parameter at the speed of light, the 
same as that of a massless particle impinging upon the Reissner-Nordstrom black 
hole. 

Although I have no real justification of Eq. (|119|) for relativistic velocities, the 
fact that it was found to be a fairly good approximation for all velocities calculated 
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so far (i.e., at least for velocities not too ultrarelativistic) for Schwarzscliild black 
holes [|lll leads me to conjecture that it may be accurate to within several percent for 
all velocities for Reissner-Nordstrom black holes as well, assuming only that /i ^ 1. 
However, apparently only numerical calculations could confirm this. 

Following Eq. (90) of p, one might note that this crude guess, Eq. ( |119| ), gives 
a minimum critical impact parameter of 
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ignoring correction terms of order /i that for /i <^ 1 would be negligible in comparison 
with the 0(1) errors I would expect to be in Eq. (|119|) . It would be interesting to 
learn from numerical calculations how the minimum critical impact parameter (say 
divided by M to make it dimensionless) , and the velocity at which it occurs, varies 
as a function of Q^/M^ for a Reissner-Nordstrom black hole. 



9 Effects of a long string on the sun or earth 

It may be amusing to conclude this hypothetical paper with a discussion of what 
would happen if a long straight string were impinging upon the sun, earth, or moon. 

Approximate these objects as spherical, with mass M, charge Q = 0, and radius 
-R, not to be confused with the magnitude of R defined in Eq. (^). Then the 
mathematical analysis above applies so long as the string does not get inside the 
object, so the separation r between the string and the center of mass of the object, 
the magnitude of r defined in Eq. (|33D, must remain greater than R. We shall 
find that for solar system objects, this puts us in the stiff-string regime unless n is 
extremely tiny, much smaller than /i ~ 10~^ expected for cosmic strings. 

In the critical capture case, inserting the units Tu and Lu from Eqs. ( |35D and 
(|36D into Eqs. (^) and (^3]), with Lq = ^cnt'^o and with &crit(wo) being given by Eq. 
(|8lD in this stiff-string regime, gives the minimum separation r^ and the relative 
velocity Vm at this separation being given by 
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(122) 



(123) 
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This expression for Vm is Jtth/A times the speed of a test mass, with the same 
charge-to-mass ratio Q/M as the object, that falls in to the surface of the object 
from a much lower velocity far away, because the static force given by Eq. ([I5|) 
between a compact object and a string is irfi/A times that between the object and 
a copy of itself at the same separation without the string. 

One can now invert Eq. ( |122|) and replace r^ by R to get the initial velocity Vq 
needed so that the critical capture occurs with the minimum separation being given 
by the radius R of the object: 

The corresponding critical impact parameter is 

For the initial velocity to obey the criterion fo ^ (1 — Q^ / M'^Y^^ jj,'^^^ for the stiff- 
string approximation to be valid during the critical capture, one needs 

(i-|;)(f )%<.«!. im 

This is an extremely weak restriction for solar system objects, since M/R is approx- 
imately 2.12 X 10~^ for the sun, 6.96 x 10~^° for the earth, and 3.14 x 10~^^ for the 
moon. 

For the sun, with mass M ^ 1.47662504 km, charge Q ~ 0, and radius R ~ 
6.960 X 10^ km, critical capture with the minimum separation being R gives the initial 
velocity and impact parameter as vq ~ 0.827(yu/10~^)^'^ m/s ~ 2.98(/i/10~^)^'^ 
km/hr and 6crit ~ 4.61 x 10"(/i/10-6)^/^ m ^ 3.08{ij/10-^y/^ AU (astronomical 
units). The relative velocity of the sun and string when the string just grazes 
the sun is v^ ~ 547(/i/10"^)^/2 m/s ^ 1970(/i/10"^)^/2 km/hr. When the string 
is grazing the sun, it gives the sun an acceleration of a = 7r/iM/(4r^) ^ 2.15 x 
10~^(;u/10^^) m/s^, which is less than the average acceleration of the earth toward 
the sun, 5.93 x 10~^ m/s^, unless /i > 2.76 x 10~^. Long slow strings with /x ~ 10^^ 
or less apparently don't give much danger of pulling the sun away from the earth, 
but they could change the orbital characteristics of the earth by a few percent, with 
possibly disastrous climatic effects. 

For the earth, with M ^ 0.444 cm, Q ^ 0, and R ~ 6371 km, one gets Vq ~ 3.65 x 
10-^{fi/m-^y^'^ m/s ^ 13.1(/i/10-^)^/^ cm/hr, 6crit ~ 1-73 x lO^^^fi/lQ-^y/^ m ^ 
11.6(/i/10-*^)i/^ AU, and v^ ^ 9.91(/i/10~6)i/2 j^/g ^ 35.7(^/10-6)1/2 km/hr, if the 
earth were isolated so that we could neglect the larger effect of the sun on the earth 
during the interaction of the string with the earth. When the string is grazing the 
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earth, it gives the earth an acceleration of a ~ 7.71 x 10~^(/i/10^^) m/s^, which is 
less than the average acceleration of the earth toward the sun if yU < 7.68 x 10""'. 
Thus a string with roughly this small a tension or smaller could not kidnap the 
earth from the sun. The average acceleration of the earth toward the moon is about 
3.32 X 10~^ m/s^, which is greater than the acceleration of the earth toward a string 
grazing it if /i < 4.30 x 10~^, so a cosmic string with fi ~ 10~^ would be unable to 
pull the earth away from either the sun or the moon (barring contrived resonance 
effects), though there is not a large margin of safety for pulling the earth away from 
the moon. 

For the moon, with M ^ 0.0546 mm, Q ^ 0, and R ~ 1738 km, one gets Vq ^ 
7.58xl0-^(/i/10-^)^/4m/s ^ 2.73(/x/10-6)^/^mm/hr, bent ~ 5.34:XlO^'^{fi/10-^y/*m 
^ 35.7(/x/10-'')i/^ AU, and w^ ^ 2.10(^^/10-^)^/2 ^/g ^ 7.58{fi/10-^y/^ km/hr, 
assuming the moon were isolated. At the moon's surface, a string would give it 
an acceleration of a ~ 1.27 x 10~^{fi/10~^) m/s^, which is less than the average 
acceleration of the moon toward the sun if /i < 4.66 x 10^^, and is less than the 
average acceleration of the moon toward the earth, 2.70 x 10"^ m/s^, if /i < 2.12 x 
10^^. Thus the moon is in even less danger of being pulled away from the sun than 
the earth is, and even if the earth-moon system were isolated from the sun, a string 
with reasonably small tension could not remove the moon from the earth. 

Another effect one can estimate is that of a string on the tides. If the earth 
had a string grazing it or running through it, presumably the curvature of space- 
time in and near the earth would be altered by a fraction of order /i. The Rie- 
mann curvature tensor at the surface of the earth has the orthonormal component 
i?%^ = 2M/R^ ^ 3.08 X 10-6 s-2 ^ 40.0 hr-^ which is about 1.79 x 10^ times the 
corresponding curvature contribution from the moon, and about 3.89 x 10^ times the 
corresponding curvature contribution from the sun. Therefore, if a cosmic string ran 
through the earth, the perturbation of the earth's tidal force would be of the order 
of 18(/i/10-6) times the tidal force of the moon. If the earth's axis were not aligned 
with the direction of the string, as the string cuts a swath through the rotating 
earth, the resulting perturbation in the earth's tidal (curvature) gravitational field 
might produce Bay-of-Fundy size tides at many coastal cities. Bids are now being 
accepted for the rights to a disaster movie. 

Acknowledgments 

I thank Jean-Pierre De Villiers and Valeri Frolov for getting me interested in this 
subject, and I thank them and Patrick Brady, Gary Horowitz, and Brandon Carter 
for further discussions. This work was supported in part by the Natural Sciences 
and Engineering Research Council of Canada. 



31 



References 

[1] J.-P. De Villiers and V. Frolov, Int. J. Mod. Phys. D7, 957-67 (1998), |r 



qc/9711045 



[2] J.-P. De Villiers and V. Frolov, Phys. Rev. D58, 105018 (1998), br-qc/9804087 . 
[3] D. N. Page, Phys. Rev. D58, 105026 (1998), |gr-qc/980408"8t 



[4] J.-P. De Villiers and V. Frolov, "Gravitational Scattering of Cosmic Strings 
by Nonrotating Black Holes," University of Alberta preprint Thy 08-98 (1998), 
|gr-qc/98T20T6| . 



[5] Y. Nambu, in Proceedings of the International Conference on Symmetries & 
Quark Models, Wayne State University (1969); Lectures at the Copenhagen 
Summer Symposium (1970). 

[6] T. Goto, Prog. Theor. Phys. 46, 1560 (1971). 

[7] B. Carter, Phys. Lett. B224, 61, (1989); J. Geom. Phys. 8, 53 (1992); Class. 
Quant. Grav. 9, 19 (1992); Phys. Rev. DIO, 4835 (1993); in Formation and 
Interactions of Topological Defects, e ds. R. Brandenb erger and A.-C. Davis 
(Plenum, New York, 1995), p. 304, |hep-th/9611054| ; "Brane Dynamics for 



Treatment of Cosmic Strings and Vortons," Tlaxcala lecture notes, 2nd Mex- 
ican School on Gravitation and Mathematical Phy sics, 1-7 Dec. 19 96, eds A. 
Garcia, C. Lammerzahl, A. Macias, and D. Nunez, |hep-th/ 9705 172 . 



[8] A. Vilenkin, Phys. Rev. D23, 852 (1981). 

[9] A. G. Smith, Tufts University Report TUTP-86-11 (1986), published in G. 
W. Gibbons, S. W. Hawking, and T. Vachaspati, eds.. The Formation and 
Evolution of Cosmic Strings (Cambridge University Press, Cambridge, 1990). 

[10] V. P. Frolov, V. D. Skarzhinsky, A. I. Zelnikov, and O. Heinrich, Phys. Lett. 
B224, 255 (1989); B. Carter and V. P. Frolov, Class. Quant. Grav. 6, 569 
(1989). 

[11] J. P. De Villiers, private communication. 

[12] B. Carter, private communication. 



32 



